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Abstract. We continue the study of constructing invariant Laplacians on Julia sets, and study- 
ing properties of their spectra. In this paper we focus on two types of examples: 1) Julia sets of 
cubic polynomials z 3 + c with a single critical point; 2) formal matings of quadratic Julia sets. 
The general scheme introduced in earlier papers in this series involves realizing the Julia set as 
a circle with identifications, and attempting to obtain the Laplacian as a renormalized limit of 
graph Laplacians on graphs derived form the circle with identifications model. In the case of cu- 
bic Julia sets the details follows the pattern established for quadratic Julia sets, but for matings 
the details are quite challenging , and we have only been completely successful for one example. 
Once we have constructed the Laplacian, we are able to use numerical methods to approximate 
the eigenvalues and cigenfunctions. One striking observation from the data is that for the cubic 
Julia sets the multiplicities of all eigenspaces (except for the trivial eigenspace of constants) are 
even numbers. Nothing like this is valid for the quadratic julia sets studied earlier. We are able 
to explain this, based on the fact that three is an odd number, and more precisely because the 
dihcdral-3 symmetry group has only two distinct one-dimensional irreducible representations. 
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1. Introduction 

This is the third of a sequence of papers (following [FS] and [ADS]) exploring the definitions and 
properties of Laplacians on Julia sets. One theme is the extension of ideas pioneered by Jun Kigami 
in the context of post-critically finite (PCF) self-similar fractals ([Ki], [S]). The Julia sets do not 
have the strictly self-similar structure, but they do share the finite ramifications properties of PCF 
fractals. Kigami's approach is to approximate the fractal by a sequence of graphs, take a sequence 
of graph energies on these graphs and then an appropriately renormalized limit to obtain an energy 
on the fractal. Combining the energy and a measure on the fractal via a weak formulation yields 
a Laplacian, which can then be described as a suitably renormalized limit of graph Laplacians, 
echoing the basic idea of calculus that derivatives are limits of difference quotients. The second 
theme is that Julia sets may be realized as circles with identifications. Indeed, the external ray 
parametrization gives a map from the circle onto the Julia set, and the dynamics on the Julia set is 
conjugate to a simple map on the circle. Thus we can hope to take simple structures on the circle 
and transport them to the Julia set. For example, the standard Lebesgue measure on the circle may 
be pushed forward to the Julia set, and this easily provides the measure we will use. The situation 
with the energy is more challenging. 

We may informally refer to this approach as the method of Peano curves. Indeed, the classical 
Peano curves also give a realization of domains in the plane as circles with identifications. Of 
course, no one in their right mind would suggest that the usual planar Laplacian be constructed 
by pushing forward structures on the circle to planar regions via a Peano curve. Nevertheless, our 
approach might be successful in the context of other fractals. 

In the first paper [FS] we studied a simple family of quadratic Julia sets, including the familiar 
Basilica and Douady Rabbit (the Basilica had been treated earlier in [RT]). In this family there 
are only a countable number of identifications on the circle, and they arise in a coherent iterative 
scheme that gives rise to a family of graphs whose vertices are identified points and whose edges are 
arcs of the circle. We used a rather ad hoc method to renormalize the graph energies in order to 
obtain an energy in the limit that is invariant under the doubling map on the circle (conjugate to 
the dynamics on the Julia set). In [ADS] the ad hoc procedure was replaced by a more systematic 
approach and this approach was used to study a few more complicated examples. 

We briefly outline the method. We start with a finite set of identified points on the circle, 
JW. To define an energy on the resulting graph we need to assign conductances (reciprocals of 
resistances) to the edges. The simplest choice would be to take the resistance to be equal to the 
length of the edge, but this is not satisfactory. Instead, we divide the edges up into different types, 
and multiply the reciprocal of the length by a factor b n that depends on the type to obtain the 
desired conductance. We then use the dynamics (in the case of quadratic Julia sets this is just the 
doubling map on the circle) to build a sequence of graphs, and energies that depend on 

the choice of weight factors {b n }. The key condition we want is that each level energy restricts to a 
multiple of the previous level energy, so E^ m \u) = r£^ m ~^(u), when u is the harmonics extension 
of u (the one that minimizes the E^ energy). This leads to a nonlinear eigenvalue problem for 
the {bn}. The key technical problem is to solve this eigenvalue problem. Note that this is exactly 
the problem encountered in the case of PCF fractals (see [P] for recent work on that case). 

In this paper we study a family of Julia sets for higher order polynomials, but with the condition 
that the polynomial has a single critical points, with canonical form z p + c. We are able to generalize 
the construction in [FS] to this case, and in particular we study in detail some examples when p = 3. 
Although the construction itself follows the same outline, the structure of the spectrum of the 
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Laplacian is quite different. One striking feature is that all eigenspaces (beyond the trivial space of 
constant functions) have even multiplicity. We will give an explanation for this phenomenon based 
on the fact that 3 is an odd number, and more precisely the fact that the dihedral-p group for p 
odd has only two distinct irreducible representation of dimension one (when p is even there are four 
such representations) . 

The other example we study is the formal mating of two quadratic Julia sets, specially a 
dendrite studied in [ADS] and the anti-Rabbit. We are able to construct the energy explicitly, and 
we examine in detail the structure on the spectrum of the Laplacian, which is more closely related 
to the spectrum for the Rabbit than for the Dendrite. We also briefly discuss two other simple 
mating for which we were unable to carry out the method outlined above. This strongly suggests 
that new ideas will be required to understand these example and other matings. 

In section 2 we describe the construction of invariant energies and Laplacians for families of Julia 
sets for p— th order polynomials with a single critical point. In section 3 we discuss matings of Julia 
sets. In section 4 we present numerical data for eigenvalues and eigenfunctions of the Laplacian for 
these examples. More data may be found on the website [ |http: / /www.math.c ornell.edu/^etotaOOl / 
]. In section 5 we discuss the structure of the various spectra, giving rigorous explanations of 
observations that were first made from the experimental results. 
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2. A Simple Family of Julia Sets 

Fix an integer p > 2, and consider Julia sets of polynomials z p + c. These polynomials have 
a single critical points, so their Julia sets are closely related to Julia sets of quadratic polynomials 
(the case p = 2). Typically the Julia sets may be realized as circles with identifications via exterior 
ray parameterizations. Rather than use the parameter c we use a parameter 9 in the circle R/Z. 
The correspondence c i->- 9 is not one-to-one, but different choices of c corresponding to the same 9 
yield topologically equivalent Julia sets. 

Assume for simplicity that 9 is rational. In fact, we will mainly be dealing with the choice of 
9 = p ( p fc_ 1 ) for k > 2. Subdivide the circle into the p arcs A n = (9 + ^,9 + ^^]. For each t in the 

circle define its kneading sequence (n , rii, ...) by the conditions (jpi>t mod 1) e A nj , j = 0,1,..., 
and identify points t, t' if and only if they have the same kneading sequence. Only rational points 
get identified, so there are only a countable number of identifications, and this gives a description 
of the Julia set as a limit of the circle with a finite number of identifications. 

For the special case 9 = p ^ fc 1 _ 1 ^ ) we observe that the k points -^trz\-> n — 0,1, ...,fc — 1 are 
all identified (kneading sequence (0,0,0,...), and these are the only identifications among points 
of the form ^zf- Among points of the form p ( p k_ 1 - ) we will identify { p p ^{ p _^ } n =i,...,k for 
each I = 0, ...,p — 1 because they have the kneading sequence (£, 0, 0, ...). This produces the first 
approximating graph with p vertices. The edges are the pairs of adjacent vertices on the 

circle. Note that p(p — 1) of the edges are loops, connecting identified points. Although the loops 
do not contribute to the energy, it is convenient to keep them in order to describe the inductive 
construction of the graph. Note that the lengths of the edges are ^^Zi] for n = 0, ft — 1. We call 

them edges of type n. For n > 1 these are loops, while for n = these are edges joining 2L±!m — - 



p(p k -i) 
and 

Inductively, we define the graph by retaining all the identifications in X^ n ' for m! — 

1, 2, m — 1 and identifying { f ' p t^ k }n=i,...,k for < I < p m — 1. These points have kneading 
sequence (no, ni, n m _i, 0, 0, ...) where I = Y^ijLo n jP' ■ Moreover, we can give a subdivision 
rule that describes the edges in X^ in terms of the edges in X^-" 1 ^ 1 ^. The edges of type n in 
Ar( m_1 ) have length pS=^wzr\^ ■ When n < k — 1 they remain edges in X^ of type n + 1. When 
n = k — 1 the edges subdivides into p—1 loops of type n for each n > 1 and p edges of type 0, 
where the sequence of types is 0, 1, k — 1, 0, 1, k — 1, 0, 1, k — 1, 0. Note that the lengths add 
up correctly: 

/ p(p-l) | p 2 (p- 1) | P k ^{p-l) \ (P- 1) = p k -\p- 1) 

Vp m (p fe - 1) p m (p fe -l) "' p m (p k -l)J P p m (p k -l) p m - 1 (p k -l) 

Figures 2.1 illustrates the identifications for the cubic Basilica (p = 3, k = 2) and cubic Rabbit 
(p = 3, k = 3). Figures 2.2 shows the corresponding Julia sets. 

We may now describe the energies on the graphs X^ that we will use to approximate the 
energy on the Julia set. Let u be a function on X^ m \ or equivalently a function on the circle that 
respects identifications, i.e, u(x) — u(y) if x and y are identified. Then, let 



(2.1) 
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Figure 2.1. Cubic Basilica and Cubic Rabbit Identifications (level 2) 




Figure 2.2. Cubic Basilica and Cubic Rabbit 




summed over all type n intervals [a;, y], and 

k-i 

(2.2) Ei m \u) = J2f>nE^(u) 

n=0 

for positive weights {b n }. Note that loops do not contribute to these energies, but we leave them 
in the sum (2.1) for simplicity of notation. Also \x — y\ = ^^fern ■ If u 1S defined on X^" 1 ^ 1 ^, we 

consider all possible extension to and let u be the one that minimizes E^K Call this the 

harmonic extension. It might seem that the harmonic extension depends on the choice of weights, 
but in fact it does not. In fact is is rather easy to describe the harmonic extension and to relate 
E l m) ((u)) to ^ m_1) (u). All the new vertices in X^ lie in the type k - 1 edges of X^ 1 ), and 
each such edge contains p — 1 sets of k identified points. We have to chose the new values of u on 
these sets of identified points to interpolate linearly between the values of « on the endpoints of the 
original edge in order to minimize the contribution to Eq from this set of edges. In other words, 
if the original edge in X^™ -1 ) was [x,y] with difference u(y) — u(x), we end up with differences 
h{u(y) — u(x)) on each of the p edges of type 0. (The other edges are loops.) So the edge that 
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contributed 



\u(y) - u(x)\ 2 

( P k - m (p-i) \ 
[ pni ) 



to E 



(m-l) 

fe-1 



now contributes 



P 



| u(v)- u ( x ) |2 

/ p- m (p-l) X 



to £^ m) , a magnification factor of p fc 1 . There are no other type edges in X^ and all the type 
n edges for n > 1 are just the type n — 1 edges in J^™ -1 ). Thus, 

(2.3) ^(«) = ^M ri ^(™- 1 )(n) 



n=0 



for the matrix 



(2.4) 



M 



/0 ••■ 
1 '• 

\o •• 

By choosing {&„} to be a left eigenvector for M, 



p fe - x \ 

: 

1 o y 



fc-i 



(2.5) 



n'=0 



(m-l). 



we obtain, 

(2.6) Ei m \u) = rEr- 1 '(u) 

Since M k = p k ^ 1 I, it is easy to see that r — p~*~ and b = (1, ri , , r~*~) give the unique 
positive solution to (2.5). So, we make this choice and renormalize by 



(2.7) 
so that 
(2.8) 
Then, 
(2.9) 



£ {m \u) 



l E, 



( m ), 



u) 



£{u) = lim £ {m \u) 



is well-defined on [0, oo] with £ = if and only if u is constant. We say u € dom{£) if £{u) < oo. 

Next we consider the behavior of the energy under the action of the map u — > u o P where 
P(x) = px mod 1. Note that P is conjugate to the polynomial z' p + c on the actual Julia set. 
Suppose E is an edge of type n in X^" 1 ^ 1 ^. Then p-E + | for £ = 0, 1, 1 are edges in ( m ) also 

of type n. Since the length of the edges are reduced by the factor p, it follows that the contribution 
from the p edges in X^ to £?n"^(uoP) is exactly p 2 times the contribution from E to S^ 1 " ^(«). 
Thus, 



(2.10) 



Et ) {uoP)=p 2 Et- 1) {u) 
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Since (2.10) is independent of n we obtain 

(2.11) E { b m \uoP)=p 2 E i b m - 1) (u) 



hence, 

2 

(2.12) £{uoP) = P -£(u)=p h ^£(u) 

r 

To define a Laplacian we need a measure on the Julia set, which we take to be the Lebesgue 
measure, /i, on the circle, since the identical points form a set of measure zero. This is called the 
equilibrium measure because 

(2.13) J foPdfx = J fdn 
We define Au via the identity 

(2.14) £(u, v) = - J (Au)vdn 

for all v € dom(£), where the bilinear form £(u, v) is obtained from the quadratic form £(u) = 
£{u,u) and the polarization identity. It follows from (2.12) and (2.13) that 

(2.15) A(uoP) =pT(A u ) p 



and this is the invariance condition of A under the action of P. 

We can also represent the Laplacian as a limit of renormalized graph Laplacians A m on the 
graphs X( m \ Note that there are exactly p m distinct vertices in and they are equally dis- 

tributed, so the right side of (2.14) is approximated by — p~ m ^ l£X(m) (Au(x))v(x). For a fixed 
, we choose v = 5 X (in other words v(x') — 5 XX >) and define 

(2.16) A m u{x)^-p m £^ m \u,5 x ) 



We can write this more explicitly by noting that each x has 2k edges (some may be counted 
twice), with two of each n type. Let e(x) denote this set of edges, and for each edge let y denote 
the other endpoint (note we may have y = x if the edge is a loop), and let n denote the type. Then 

(2.17) A m u(x) = ^ c n(u(y) - u(x)) 



with the conductances c n given by 
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3. Formal Matings of Julia Sets 

Suppose J' and J" are Julia sets for the polynomials P(z) = z p + d and P'(z) — z p + c" 
with the same degree p. If we realize J' and J" as circles with identifications corresponding to 
the parameters & and 9" , then we can realize the formal mating J of J' and J" as a circle with 
identifications t\ = t 2 if and only if t\ and t 2 are identified in J' or 1 — 1\ and 1 — t 2 are identified in 
J" . The main problem we consider is how to construct an invariant (under t — > pt mod 1) energy 
on J if we know how to construct invariant energies on J' and J". More specifically, suppose 
we have increasing sequences X'^ and X" {m ^ of identified points in J' and J". Then we form 
X( m ) = X'( m ) U (1 — X"( m ^) and make the appropriate identifications. We will always assume that 
X'( m ) = pX'( m -^ mod 1 and X"^ = pX" {m -^ mod 1, so that X^ = pX^ m ^ mod 1. We 
make X^ into a graph by declaring edges between all neighboring points on the circle. 

Our strategy would be to define different types of edges in these graphs, then form the type n 
energy via (2.1) and the full energy via (2.2) for a suitable choice of weights. We would then want 
to define subdivision rules for how an edges of type n in X( m_1 ) breaks up into edges of different 
types in X^ m K We could then find the energy minimizing extension u and compute a matrix M for 
which (2.3) is valid. Note that we can expect that the matrix M depends on the weights {b n }, so 
that (2.5) becomes a nonlinear eigenvalue problem. Nevertheless, if we are able to find a positive 
solution, then we can define a renormalization (2.7) of the energy (2.2) such that (2.8) holds, and 
then define an energy on J via (2.9). 

It appears that this strategy is successful in only a few special cases. We will discuss one 
example where it works, and then very briefly two examples where we were unable to make it work. 
In all these examples p = 2, so we are considering formal matings of quadratic Julia sets. From the 
point of view of dynamical systems, the most interesting formal matings are the ones which yield 
Julia sets of rational mappings. These are called simply matings. Our main example is not such a 
mating. See [] for details. 

Our main example has 0' = \ and 9" = |, so J' is the Dendrite discussed in detail in [ADS] 
and J" is the anti- Rabbit (the mirror image of the Rabbit p = 2, k = 3). In this example X'^ 
and 1 — X"( m ) consist of exactly the same points, only with different identifications. Thus X^ 
consists of all points of the form f^fery for n = 1, 2, 3 and < £ < 2 m - 1. 

The rabbit identification simply identifies j^l^ , ^{2^) anc ^ 7(2™ ) • ^ e wm ca ^ tnese triplets. 
The Dendrite makes many of the same identifications, but also makes some different identifications 
that result in pairs of triplets begin identified to create sextuplets. The first example occurs when 
m = 4, with the identification of the I = 1 and i = 9 triplet. Therefore, for m > 4, there will 
be identifications of pairs of triplets for all £ = 1 mod 2 3 . For m = 6 a different type of sextuple 
arises by identifying the £ = 5 and £ — 7 triplet, leading to identifications for all £ = 5 mod 2 5 
for m > 6. Note that the triplets first identified for m = 4 and m = 6 contain the dividing points 
j2 and J2 for the dendrite. More generally, for m = 2 j + 2 there will be triplets corresponding to 
£ = |(2 2j — 1) and £ = |(2 2j — 1) + 2 2j+1 that contain these dividing points and so get identified in 
a sextuplet, and for all m > 2j + 2 this will result in identifications of all triplets with £ = g(2 2j — 1) 
mod 2 2j+1 . The exact description of the pairing of such triplets is rather complicated and will not 
be given here. These give all the sextuplets, and if we denote by a m the number on level to, then 



2o>m 

2a m + 1 



if to is even, 
if to is odd. 
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which leads to, 




3 

2-4" 



for m > 2 

The graph X^ 11 * 1 consists of intervals of lengths 



7(2-) 



for j — 1,2,3, and we will call these 



intervals of type 1, 2 and 3. When we pass to the next graph X^ n+1 \ the intervals of type 1 
and 2 do not subdivide and simply become intervals of type 2 and 3, respectively. Intervals of 
type 3 subdivide into four intervals of types 1, 2, 3, 1 in that order. Some intervals will be loops 
(endpoints identified), and those do not contribute to the energy. Those that are not loops will 
always be paired, meaning that there will be two (and sometimes four) intervals with identified 
endpoints. Figure 3.1 shows the subdivision of paired intervals in X^ when increasing m to 5, 
yielding two sets of paired type 1 intervals. A more detailed analysis of identifications in subdivision 
will be given in section 5. Here the important observations is that the intervals of type 2 and 3 
created in the subdivision are loops and so do not contribute to energy, and intervals of type 1 are 
all paired. 



Figure 3.1. 



16/112 15/112 



9/112 8/112 




64/112 65/112 



71/112 72/112 



The analysis of the energy is then exactly the same as for the Rabbit. We consider the three en- 
ergies E[ m \u), E ( ™\u) and E { ™ ] {u) and observe that E ( ™\u) = 4 m_1) («), Ef '(«) - ^ m_1) (u) 



and E[ m \u) 



Thus, 



4£ ; g Tn so the matrix is 



M = 




(3.1) 



- 4"f (E[ m \u) + ^E { ™\u) + 4l^ m) (n)) 



Perhaps the simplest nontrivial matings are the Basilica/Rabbit and Rabbit/Rabbit matings. 
We studied these examples extensively searching for a way to sort the edges into a finite number of 
types, and then find subdivision rules for passing from one graph to the next. What we found was 
that each time we increased the level m, it was necessary to add new edge types. It seems unlikely 
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that there would be only a finite number of edge types that would work for all m, although there 
is the remote possibility that we stopped the search too soon. We then investigated constructing 
energies using a small number of edge types and weights chosen through trial and error, attempting 
to obtain approximate Laplacians that had some rough consistency from level to level. The results 
were not very satisfactory. Some data from these attempts may be found on the website. 



LAPLACIANS ON JULIA SETS III 



11 



4. Numerical Data 

In this section we present numerical data on the examples above. In each case, to obtain 
the data, we approximated the Julia set, or formal mating of Julia sets, by a graph. We then 
computed the graph Laplacian on the approximating graph, with the adjustment to the standard 
graph Laplacian coming from the weighting of the edges of the graph, and the measure placed 
on the graph. Using the built in eigenvalue and eigenvector functionalities in Mathematica and 
MATLAB, we produced approximations of the eigenvalues and cigenfunctions of the Laplacian on 
the set. 

The data presented comes from the highest level approximation computed. We present a 
selection of eigenfunctions, the first 150 eigenvalues, the eigenvalue counting function and the Weyl 
ratio. 

The first example, the cubic Basilica, corresponds to p — 3 and k = 2, so 9 = the first Julia 
set shown in figures 2.1 and 2.2. Figure 4.1 shows the first eight non-constant eigenfunctions, figure 
4.2 shows the eigenvalue counting function, N(t) = #{Aj : Xj < t}, and figure 4.3 shows the Weyl 
ratio W(t) = N(t)/t a for the appropriate power a, which in this case is a = |. Table 4.1 shows 
the first 150 eigenvalues. The second example, the cubic Rabbit, corresponds to p — 3 and k = 3, 
so 9 = jg, the second Julia set show in figures 2.1 and 2.2. The corresponding data here is shown 
in figures 4.4-6 (with a = |) and table 4.2. 

The last example is the dendrite-antiRabbit formal mating discussed in section 3. Figure 4.7 
shows the first eight eigenfunctions and figure 4.8 shows a multiplicity 5 eigenspace. Figure 4.9 
shows the eigenvalue counting function and figure 4.10 shows the the Weyl ratio with the value 
a = X determined experimentally. Table 4.3 shows the first 150 eigenvalues. 



Cubic Basilica Data 

Figure 4.1 Eigenfunctions Corresponding to First 8 nonzero Eigenvalues 




Figure 4.2 Eigenvalue Counting Function 



7 I 1 1 1 1 1 1 1 r 




n I i i i i i i i I i 

01234567B9 10 



Figure 4.3 Weyl Ratio, a = 2/3 
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Table 4.1 First 150 Eigenvalues 
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Figure 4.5 Eigenvalue Counting Function 



Table 4.2 First 150 Eigenvalues 




Figure 4.6 Weyl Ratio, a = 3/4 
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Figure 4.7 Eigenf unctions Corresponding to First 8 nonzero Eigenvalues 




Figure 4.8 5-Dimensional Eigenspace Corresponding to the Eigenvalue 893.5 



■ 

OS 

-LJ 



If, 



I ■ 

I' 

OS 

~ 

as 

-IS 



_1 i i I i i i I i i i I i i i I i i i L_ 



I ■ 



-OS 
-t.C 
-L.5 



I D 



_1 i i I i i i I i i i I i i i I i i i L 



02 OA C£ OK 



-in 



I 3 - 
lO- 
OJ - 



I, 



-0J 
-IB 



Figure 4.9 Eigenvalue Counting Function 



Table 4.3 First 150 Eigenvalues 



LDO : 
■■ - 

I'- 



r 



500 in 



3 x L!l1 sc Ltf 



Figure 4.10 Weyl Ratio, a = .7 
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5. Properties of the Spectrum 

We consider first the simple Julia sets discussed in section 2 with parameters (p, k). Note that 
(2.15) implies that if u is an eigenfunction of —A with eigenvalue A, then u o P is an eigenfunction 
with eigenvalue p^ A. In particular, the spectrum is preserved under multiplication by pT. It 

fc-Ll 

therefore makes sense to split the nonzero eigenvalues into the primitive ones, not of the form p~*~ A' 
for some other eigenvalue A', and the derived ones that have this form. Moreover, the derived 
eigenvalues have anorder j, the unique value for which the eigenvalue for which the eigenvalue 
pp~^\' for A' a primitive eigenvalue. Note that primitive eigenvalues will have order zero under 
this definition. Let to(A) denote the multiplicity of the A-eigenspace 

Conjecture 5.1. Let p = 3. Then to(A) only depends on the order of X, m(A) = mj where 
too = 2 and the recursion relations 

(3(mj_i - 1) + 1 if j = 1,2, ...,k- 1 mod fc, 
5.1 rrij = < , 

v ; 3 \3(mj_i - 2) +2 ifj = mod k. 

hold. In particular, all multiplicities are even (for A ^ 0). 

We will sketch an argument for why we believe the conjecture is valid. It is unlikely that this 
can be made into a proof since it uses the "no coincidences" principle: unrelated eigenfunctions will 
have different eigenvalues. (A similar principle may be invoked in describing the spectrum of the 
ordinary Laplacian on the unit disc, where it amounts to saying that Bessel functions of distinct 
orders do not have common zeros.) The occurrence of coincidences would increase the multiplicities 
(but they would still be even). Because of the spectral clustering phenomena discussed below, we 
expect that there will be many near coincidences, so numerical approximations will not be effective 
in ruling out coincidences. 

We note that the conjecture gives a description of multiplicities that differs substantially from 
the examples with p = 2 discussed in [FS]. We expect that suitable modifications of the conjecture 
will be valid for odd values of p. 

The key observation is that there is a symmetry group of dihcdral-3 type acting on these Julia 
sets. It is clear that the rotation t — > t + | is a symmetry, and this generates a Z 3 symmetry 
group. When k = 2 it is also clear that the reflections of the circle about the diameters (|, |), 
(t5> 15 )' (t5' 15 ) are symmetries of J . For k > 3 there are "reflection" symmetries that are not 
true reflections of the circle, but still fix one of the points in and permute the other two. For 
example, this map will actually be the identity on the whole arc containing the identified points on 
the circle corresponding to the fixed points in in other words on the set of "ears" of J that 

attach to the point. A similar idea was discussed for p = 2 in [FS]. 

Symmetries preserve the eigenspaces, so each eigenspace splits into a direct sum of spaces that 
transform according to the irreducible representations of the symmetry group. In this case there 
are just three such representations: two one-dimensional representations, the trivial and alternating 
representations, and one two-dimensional representation. Note that both one-dimensional repre- 
sentations act like the identity on the rotational Z3 subgroup. 

Lemma 5.2. An eigenfunction u is of the form v o P for another eigenfunction v if and only if 



(5.2) 



«(*+ 3) = «(*) 
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,4s a consequence, a primitive eigenspace cannot contain any components corresponding to one- 
dimensional representations. 

Proof. If (5.2) holds, we may define v(t) = Then u(t) — v(3t) and v is an eigenfunction. 

Conversely, if u(t) — v(3t), then clearly (5.2) holds. □ 

This shows that primitive eigenspaces have even multiplicity, and the no coincidences principle 
would make mo = 2. However, we can arrive at another important insight with the observation that 
the two dimensional representation contains one dimensional subspaces that are skew-symmetric 
under one of the reflections. Thus, the eigenspace contains a function u that satisfies u o R = —u 
for the reflection that fixes a point x in Now x splits J into a large arc and a small arc, 

and because of the skew-symmetry we can write u as a sum of eigenfunctions supported on each of 
them. However, we may argue that the one supported on the small arc must be zero, for if not we 
could rotate it around to create a nonzero eigenfunction satisfying (5.2), contradicting the fact that 
we have a primitive eigenspace. Therefore we conclude that the primitive eigenspace has a function 
supported on the large arc and vanishing on the small arc. Taking two points in gives us a basis 
(of course the functions associated with all three points are linearly dependent). We already observe 
a contrast with the p = 2 case discussed in [FS] where essentially every eigenfunction vanishes on 
an open set, since if we take a sum of two such functions the result will not have such a vanishing 
property. 

Next we describe how to pass from a primitive eigenspace to the derived eigenspace of order 
1. Let u be the eigenfunction that vanishes on [ 3 fe_ 1 > 3fc_i ]- Then v{x) = u(3x) vanishes on 

[ 3 (3fc-i) ' 3(3* -1) 1' [ 3 (3fc-i) ' ^t^-i)" 1 ] and [ 3 2 (3*-i) ' 2 ' 3 3(t*-"i)~ 2 ]- Notc that v is supported on the 
union of three intervals / = [§k^, §e3t]j I + \ an d I + §, and these intervals are separated in J. 
Therefore we can split v into a sum of three eigenfunctions supported in each of these intervals. On 
the other hand v(x + |) = u(3x + |) is supported on I + |, I + | and I + | and these intervals 
are connected in J . Thus we arrive at the conclusion that mi — 4 with a basis of eigenfunctions 
vanishing on certain open intervals. 

The process of passing from derived eigenspaces of order j — 1 to derived eigenspaces of order 
j is similar. If (ui, ...,u mj l ) is our special basis for the order j — 1 eigenspace, we consider the 
functions v m (x) — u m (3x). Depending on the residue of j mod k, all but one or all but two of 
them will be supported on separated intervals, and may be split into three linearly independent 
eigenfunctions. This gives the relation (5.1), and a special basis. 

We can make another conjecture concerning the location in the spectrum of derived eigenval- 
ues. Let A = denote the lowest eigenvalue, and then label consecutive eigenvalues Ai, A2, A3, ... 
repeated according to multiplicity. For any eigenvalue A denote by #(A) the largest n such 
that A„ = A. We know that the eigenvalue counting function N(t) has growth rate t a for 

k 

a = '"to^* 1 = £qpj and the Wcyl ration tends to a multiplicatively periodic function of 

fc-f- 1 1 k I 1 

period p~^ . That means #(p J_ ^^) ~P : '#(^)- 
Conjecture 5.3. Letp = 3. Then 

(5.3) #(3^A)=3^#(A). 

This conjecture is based solely on experimental evidence. It is not valid for p — 2. A similar 
conjecture was described in [ADS] for the Dendrite example, but not for the other examples pre- 
sented there. A similar result was proved for the Sicrpinski gasket in [S]. The experimental evidence 
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points to the existence of spectral gaps between A^-i and A3-, although it is not clear what the 
lower bound of the ratio A ^~ 1 might be. Also, we see many eigenvalues that are distinct rela- 
tively but close together, suggesting that spectral clustering (arbitrarily many distinct eigenvalues 
in arbitrarily small intervals) may occur. 

Next we discuss the spectrum for the Dendrite- Anti Rabbit mating. For this we need a precise 
definition of the two reflection symmetries. For the rabbit there are reflections we called horizontal 
p H and vertical py in [FS], so we will keep the same terminology here, since the reflection on 
the mating are essentially the same. (Note that these are not symmetries of the Dendrite.) The 
horizontal reflection is the map t — > t + \ on the intervals [^, |] and [yy, yy] and maps [yy, ^] to 
itself, and also [|, yy] to itself, while the vertical reflection is the identity on the intervals [|, |] and 
[yy, yy] and interchanges the intervals [yy, ^] and [|, yy]. 



Figure 5.1. 






To describe those reflections on the remaining intervals [yy, |] and [|, yy] we observe that we 
may define an "orientation reversing" map between two intervals in the graph of type 3 

iteratively as follows. Say the intervals are [f^fSj, 7(2™) ] an d [fp^)' 7(2™ ) ] ' ano - we wairt to map 

T(0j t0 T(0) and Y(2^) to Ip^J ■ We subdivide the intervals in X (m+1) and map the outer 
subintervals to each other in orientation reversing fashion, and the inner subintervals are mapped 
to each other by translation. See Figure 5.2. In particular, pn is the identity on [yy^, yy|] and maps 
[1T2' TT2] to Ink' Til] reversm g orientation (and similarly on [f , yy]), and py is the translation 
t — > t + \ mapping [yy^ , yyg] to [yy^, yyj] (and backward), while it is orientation reversing between 
[yfg, yfj] and [yy^, yy|] and similarly between [y\|, y^] and [yy|, f^]. 

These two reflections generate a symmetry group isomorphic to Z 2 x Z 2 , the composition in 
either order being the translation t — > t + |. Primitive eigenfunctions must be skew-symmetric 
with respect to this translation, so they fall into two types: horizontal eigenfunctions are skew- 
symmetric with respect to pn and symmetric with respect to py, while vertical eigenfunctions have 
the reverse symmetry. Note that horizontal eigenfunctions vanish on [yyjj j^] and [jf2, j^} because 
Ph is the identity on those intervals. It then follows that eigenfunctions derived from horizontal 
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Figure 5.2. 
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cigcnfunctions will vanish on prescribed intervals, but will generate one-dimensional eigenspaces 
(assuming no coincidences). 



Similar reasoning shows that vertical eigenfunctions are supported in 



tually we can say more. 

Theorem 5.4. Primitive vertical eigenfunctions are supported in 
t then u(pv{t)) = —u(t) imples u{t) = 0. 
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Proof. If pv(t) ~ t then u(pv{t)) = —u(i) imples u(t) = 0. From this we infer that u is 
supported on[^,^]U[|,^], and it also vanishes at the sextuplet n2^iT2^Ti^~ii2~~ii2~~ 
It also follows from the skew-symmetry that we may cut the eigcnfunciton at these points, 



71 
112 

to obtain three distinct eigenfunctions supported on 
9 1L] u r_g5_ _67i and on \JJ_ A5_] ij r_67_ jn_i 
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We need to show that the first two vanish 



identically, and for this we need to use the fact that we are assuming the eigenvalue is primitive. 
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is also an eigcnfunction, contradicting the assumption that the eigenvalue is primitive. Similarly, if 
u is supported in the first union of four intervals, 



v(x) 



u{\x) 
-u{\x- 



16 ) if a; € [-j-j, j] 
otherwise. 



is also an eigenfunction. 



□ 



U 



We can now describe the derived eigenspaces and their multiplicities. If u(t) is a primitive 

vertical eigenfunction, then u(2t) is a derived eigcnfunction supported on , ? / 5 ^ 1U[ 

123 127 1 U [yM), 7py]- This support is not connected 
15 
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however, and splits into the connected 
1 1 ' 1 By restricting u{2t) to each of these 
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sets we obtain a 2-dimensional eigenspace. By iterating this argument we obtain a splitting of 
u(2 k t) into 2 k linearly independent eigenfunctions. But this is not the whole story. There are two 
other types of eigenfunctions that arise. 

To describe the first type consider the case k = 2. Each of the four pieces of u(4t) is supported 
in a pair of intervals of length ypiy, and the eight total intervals are obtained from [^pey , yp^y] by 
rotation by multiples of |. However, there is another pair of intervals [ 7 / 3 2 9 6) , ypgy] and [fpty, fpFy] 
that are paired and support another eigenfunction obtained from pieces of u{4£) that are rotated 
by multiples of ^ . Thus we get a multiplicity 5 eigenspace. When k = 3 we double the multiplicity 
to 10 by splitting each of the 5 eigenfunctions from k = 2. For k — 4 we split the 10 previous 
eigenfunctions, add one more by paired intervals rotated by multiples of but also there is yet 
another type of eigenfunction supported on a set of four intervals with endpoints belonging to 



two sextuplcts of identified points 



36 
7(2 6 ) ' 



37 
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This is 



7(2 6 )J w L7(2 6 )' 7(2 6 )J w L7(2 6 )' 7(2 6 )J w L7(2«S)> 7(2<3)J 

illustrated in Figure 5.3. This is the j = 2 case of a more general construction that yields an 
eigenfunction supported on 2-? intervals with endpoints belonging to 2- 7 ~ 1 sextuplcts of identified 
points. These eigenfunctions arise at level k for k > 4 even with j < | . 



Figure 5.3. 
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In general, when k is odd we just split the eigenfunctions from level k — 1. When k is even (and 
k > 4) we also add a total of | eigenfunctions are described above. If bk denotes the multiplicities 
then 



26 fc _i 
26 fc _i + 



if k is odd, 

if k is even, and k > 4, 



This implies b 2 k = 46 2 fc-2 + ^> which is easily solved to obtain 6 2 fe = ip4 fe — | — |, and then 
W=2(f4*-§-§). 
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